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Abstract 
In this work, we employ renormalization group methods to study the general behavior of field theories 
, possessing anisotropic scaling in the spacetime variables. The Lorentz breaking symmetry that accompanies 

O^l ■ these models are either soft, if no higher spatial derivative is present, or it may have a more complex structure 

q; if higher spatial derivatives are also included. Both situations arc discussed in models with only scalar fields 
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and also in models with fermions as a Yukawa like model. 
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I. INTRODUCTION 

Quantum Field Theories with anisotropic spacetime scahng have been recently considered in the 

n 

literature in the search for viable frameworks for quantum gravity [1] and also for the completion 
of nonrenormalizable effective field theories [4] . The fundamental assumption of theses approaches 
is that time and space behave differently under a general scaling: x* — >■ 62;* whereas t -^ b^t. 
With the "critical exponent" z conveniently chosen, renormalizability can be achieved without 
the undesirable introduction of ghost degrees of freedom which would be present if higher time 
derivatives were also included [3|. Various investigations on the properties and applications of 
these approaches may be found in the literature [4]. However, the anisotropic scaling inevitably 
encompasses a breaking of Lorentz invariance so that a basic question here is if and in what 
circumstances this symmetry is restored. Although an intuitive argument says that the restoration 
may take place at low momenta, a general study of the renormalization group flows for scalar 
models indicates that in reality it requires a careful analysis [5| (see also a]). 

Concerning the Lorentz symmetry breaking, different situations may arise from modifications on 
the kinetic part of a given Lagrangian that keep the dependence on the time derivative unchanged: 

1. The coefficient of the spatial derivative is modified but higher derivatives are not present. In 
models for a single field this is physically innocuous since it may be adjusted to any finite value 
by a mere change of dimensional units; no breaking of Lorentz symmetry really occurs. However, 
in models with more than one field the modification may have physical implications with radiative 
corrections which propagate with^different velocities such that the breaking of Lorentz symmetry 
may increase with the energy. In [S] this possibility was shown to occur in a model of fermions and 
bosons interacting via a trilinear coupling in four dimensions. As we shall see, the same happens 
in a model of two scalar fields coupled through trilinear interaction terms in six dimensions. For 
completeness, we also discuss the case of a Yukawa like model. Whenever the Lorentz symmetry 
may be restored we say that the breaking is soft. 

2. Higher spatial derivatives are introduced, i.e., we have a truly anisotropic model. In this 
case, the coefficient of some of the higher derivatives generally changes by effect of the radiative 
corrections in such way that the breaking of Lorentz symmetry is either soft in the infrared or 
becomes stronger by decreasing the energy. It must be stressed here that by soft in this anisotropic 
situation only means that the effective coefficient of a higher spatial derivative term decreases; in 
general it can not vanish as the model may become nonrenormalizable. In this situation, we shall 
consider the possibility of the Lorentz symmetry restoration in an approximated sense. 



In this work we will pursue these studies focusing some aspects of the renormalization of 
anisotropic field theories. With the ultraviolet improved free propagators, usual Lorentz sym- 
metric models that are renormalizable in a certain dimension, will become renormalizable in a 
higher dimension. Thus with z = 2, the (p^ and 99^ models that are renormalizable in 6 and 4 
space-time dimensions will become renormalizable in 10 and 6 spatial dimensions, respectively. It 
should be noticed that although closed forms for the Feynman amplitudes are in general unfeasi- 
ble, it is possible to calculate the renormalization constants and to determine the flows of relevant 
coupling constants. This may be done by an application of the BPHZ renormalization theorem to 
dimensionally regularized integrals as it will be described shortly. As an application, we consider 
the ifi^ model in ten dimensions and with z = 2 which is asymptotically free. Due to this property, 
it is not possible to obtain reliable results in the small momenta region. Nevertheless, this model 
furnishes a simpler setting to expose our methods and discuss some properties of anisotropic field 
theories as unitarity. Afterwards, we analyze the tp'^ model in six space dimensions which presents 
an increased degree of complexity as the calculations of the renormalization constants involve a 
two-loop diagram. In spite of that, for a special configuration of the initial values of the effective 
parameters, we found an explicit solution for the renormalization group equations. 

We also consider a model of fermionic and bosonic fields interacting via a Yukawa like Lagrangian 
in six space dimension with z = 2. As we shall see, differently from the ip^ model in 10 spatial 
dimensions, this model is infrared stable at the origin and so it seems to be a good candidate for 
testing the recuperation of Lorentz symmetry at low momenta. Indeed, for the four dimensional 
model and without the higher derivative term we will demonstrate that Lorentz symmetry is 
restored for small momenta if initially it was broken by assuming different light velocities for the 
boson and fermionic components. When higher spatial derivatives are introduced, we performed 
a numerical analysis, which is required by the complexity of the renormalization group equations 
which prevents the existence closed analytical solutions. In this last situation, we found that the 
Lorentz symmetry breaking parameters also decrease in the infrared region. 

Our work is organized as follows. In section II, we consider a six dimensional model of two 
scalar fields interacting through trilinear couplings as a tool to study the soft breaking of Lorentz 
symmetry; in this preliminary study no higher derivatives are present. In this context, we also 
analyze the soft breaking of the Lorentz symmetry in a Yukawa like model. In section III, we 
present some general remarks and the formalism we shall use in this work. A subsection is devoted 
to the unitarity problem in higher spatial derivative models. Sections IV and V are dedicated 
to the analysis of renormalizable versions of the p'^ and Yukawa models which are infrared sta- 



bles. A summary and additional remarks are presented in the Conclusions. Two appendices to 
study integrals needed in our studies and provide additional details on the computation of the 
renormalization group parameters are included. 

II. SOFT BREAKING OF LORENTZ SYMMETRY 

Before considering truly anisotropic models which contains higher derivatives here we will an- 
alyze the soft breaking of Lorentz symmetry in two situations, namely, a six dimensional model 
with two scalar fields and a Yukawa like model in four dimensions. Let us begin by considering the 
purely bosonic model described by the Lagrangian density 
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The free propagators for the 99 and (j) fields derived from this Lagrangian, 
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will be represented by dashed and dotted lines, respectively. With this graphical notation, the one- 
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loop contributions for the two point vertex functions F^ and r\ are drawn in Fig. [TJ Notice that 
the corresponding analytic expressions, which are supposed to be regularized by taking the model 
to D = 6 — e dimensions, are exchanged by the replacements ((/?, Ai, A4, A3) ^ [(j), A2, A3, A4). The 
renormalization group flows of the parameters in ([1]) are flxed by the introduction of dimensionless 
coupling constants Aj — > /U2 Aj and the computation of the pole part (PP) of the relevant diagrams. 
Following standard procedures, we found that the dependence on the external momentum of 
the pole part of the diagrams in Fig. [1] are 
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whereas PP[Fi^i!\ and PP[Fi^^] are obtained from PP[Fi^;)\ and PP[Fi^p\ by replacing A4 
by A2 and Ai by A3, respectively. Using these results we may obtain the beta functions associated 
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FIG. 1: One-loop corrections to the two point vertex functions of the tp (first row) and (p (second row) fields. 



with the parameters 6^ and 6^: 
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The basic reason why these expressions do not depend neither on Ai or on A2 is that the would 
be contributions are the same as in the case where there is only one field, either (p or (f). As we 
mentioned in the Introduction, in those cases the parameter h^p or 6^ does not receive radiative 
contributions what implies that the corresponding /3 function vanishes. 

Notice that b^ = b^p ^ b is a. renormalization group fixed point where the Lorentz symmetry 
holds. Now, if b^p > b^f, then, by lowering the energy, 6^ decreases and 6^ increases and the fixed 
point is infrared stable, the two parameters approaching the fixed point value. Similarly, if 6,^ > 6,^, 
b(f, decreases and b^p increases by lowering the energy the two parameters again tending to the fixed 
point b. We conclude that in any case the Lorentz symmetry is restored in the low energy regime. 

We will also determine the pole part of the three point vertex functions. Observe that in this 
case any given diagram either has all three internal lines of the same type or has two equal lines 
but of type different from the remaining one. In the last situation the integrals to be computed 
have the general form 
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where A^ symbolically represents the factor associated with the coupling constants and either 

{A,B) = {b^,b^) or {A,B) = (b^^b^p). A straightforward calculation yields then the following pole 



part: 
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The above result allows us to calculate the beta functions associated with the various coupling 
constants. For example, we obtain 
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Let us now turn our attention to a Yukawa like model specified by 
1 52 ^2 _ 

+ i^v^r^^v. - ^y. (9) 



In this case we find 
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and 



which, analogously to the purely bosonic model, shows that b^ = b^ is also an infrared renormal- 
ization fixed point where Lorentz symmetry holds. 



III. ANISOTROPIC SCALING: GENERAL CONSIDERATIONS 

The anisotropic field theories with bosonic and ferniionic components that we will consider have 
the generic form 

^ = ^0 + ^int, (12) 

where the free part, Cq, is given by 
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for bosonic fields and 
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for fermionic fields. In these expressions each latin index runs from 1 to d, the spatial dimension of 
the model and the effective signs of the a's and /3's have to be chosen so that the energy associated 
with Cq is positive. Cint describes the interaction between these fields. Notice that z designates 
the highest degree of the spatial derivatives and we have included terms with less derivatives as 
they may be necessary in the renormalization process. As each power of xq scales as z powers of Xi, 
the effective dimension of the Lagrangian is z + d. Then, by taking a^ and (3^ to be dimensionless, 
we find that effective dimensions of ip and ip are respectively 

Dim[c^] = and Dim[V'] = — • (15) 

We will be dealing with Feynman amplitudes of the form 

L 
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where L is the number of loops, k = (/ci, . . . , /c^) and p = (pi, . . . iPn) are, respectively, the loop 
momenta and external momenta associated to a generic proper (i.e. IPI) Feynman diagram G. The 
unsubtracted amplitude Iq is a product of anisotropic propagators and monomials in the momenta 
of the lines joining at the vertices of G, 

lG{k,p,m) =WPa{k,p)\[AF{labc), (17) 

a abc 

with 

P{1) 

where / = (Iq,^) is the momentum flowing through a line of G, the bi are simple functions of the 
a's , /3's and M (bi = ai for a bosonic propagator and m = M for a fermionic one), and P{1) = i 
for a bosonic propagator and a polynomial of first degree in l^ and of z degree in 1, for a fermionic 
propagator. The degree of superficial divergence of G is fixed by power counting in which each 
power of the time like component of a vector counts as the z power of its spatial like component. 
This gives 

d{G) = id + z)L-2znB-znF + ^Da, (19) 

a 

where hb and np are the number of bosonic and fermionic internal lines of G; Da is the degree of 
the monomial Pa{k,p) assigned to the vertex Va- Now, if G has V vertices, L = ns + np — V + 1 



and therefore 

d{G) = d + z-{z- d)nB + dnp + Y^{Da -d-z). (20) 
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We also have the topological identities 

2nB + NB = ^Va and 2nF + TV^ = ^ i^f , (21) 

a a 

where Nb and Np are the number of external bosonic and fermionic lines, v^ and v^ are the 
number of bosonic and fermionic lines joining at the vertex Va- Using these relations, we obtain 

d{G) = d + z- V)\m.[ip\NB - Dim[V']iVF - ^{d + z- Dim[K]), (22) 

a 

where Dim[V^]= L'a+Dim[(^]f^+Dim['j/^]f^ is the "canonical" operator dimension of the term in 
Cint associated with Va- As usual, we classify a given vertex Va as being non-renormalizable, renor- 
malizable or super-renormalizable according it has dimension greater, equal or less than z+d. Thus, 
a purely fermionic theory with a quartic non-derivative self-interaction is renormalizable \i z = d 
and super-renormalizable \i z > d. A renormalized amplitude associated with the graph G may 
)e obtained by applying subtraction operators arranged according Zimmermann' s forest formula 
7|. For the special case in which G is primitively divergent (i.e., without divergent subgraphs), 
the integral in (J16p can be made finite by replacing I{G) by 
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where [x\ is the greatest integer less than or equal to x, pg symbolically stands for the product of 
s time-like components of an independent set of external momenta; pi denotes the i-th space-like 
momentum (with the index of the component implicit) and all derivatives are computed at zero 
external momenta. Actually, in our one-loop calculations performed in the next sections we will use 
the above result just to unveil most easily the pole part of dimensionally regularized amplitudes; 
afterwards we apply our renormalization prescription which consists in removing these pole parts, 
what is usually called MS subtraction scheme. Besides, throughout this work we will take the 
critical exponent z to be two. 

As in the usual isotropic situation, massless theories requires a special consideration. In these 
cases it is better to use modified Taylor operators so that the last subtraction is performed replacing 
m by an auxiliary mass parameter /i which plays the role of a renormalization point. 

Using this scheme, the BPHZ normal product algorithm can be extended to the present situ- 
ation. Thus, if O is a formal product of the basic fields and their derivatives, a normal product 



of degree 5, Ns[0], is defined in the usual way [8]. Notice that 5 = Dim[0] + c, where c is a 
non- negative integer and the dimension of O, Diui[0], is computed counting z for each "time" 
derivative, 1 for each spatial derivative and the dimensions of the basics fields as fixed before. As 
in the isotropic situation, these normal products satisfy a number of convenient properties which 
allows a systematic way for deriving Ward identities and computing their anomalies. In particular, 
it should be noticed that, inside a Green function, 

diNs [O] = Ns+i [diO] , doNs [O] = Ns+, [doO] . (24) 

A simple example 

As a simple example of the methods exposed in the previous section, we will examine a renor- 
malizable anisotropic (p^ model which was first considered in ^. Renormalizability requires 
Dim[99^] = d + z, which for z = 2 fixes d = 10. In this situation, the degree of superficial di- 
vergence of a proper graph G is given by 

d{G) = 12 - AN, (25) 

where A^ is the number of external lines of G. Taking this into consideration, the Lagrangian 
including counterterms to cancel the pure pole part is 
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where A is the spatial Laplacian and 
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with the constant c adjusted to eliminate all tadpoles. Thus, at one-loop order divergent graphs 
have two or three external lines and are quartically or logarithmically divergent, respectively. 

Let us analyze each case separately: 

I. One-loop correction to the two point function. The regularized amplitude is given by 
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where k"^ = k^ — b'^k.'^ and d = 10 — e. As mentioned before, because of the higher power of the 
spatial momentum in the denominators, S does not possess a closed analytic expression. However, 



to compute its pole part it is enough to calculate the action of the Taylor operator f^ on the above 
integrand. This leads to integrals which, after using the spatial rotational symmetry, have the 
general form 

f dko d<^k k^\k\y 

where for d = 10 the parameters x, y and z are such that 2x + y — Az + 12 is either equal to 0, 2 or 
4. In these cases, neglecting finite parts, J{x,y,z) is given by Eq. ()A6p in the Appendix lAl With 
the help of that result, we have: 

a. Term with four derivatives. Omitting some finite contributions, the term with four spatial 
derivatives yields 
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where den = fc^ — m? — a^(}i^)'^ + ie . Performing the integrals, we find that their pole part yields 



the result 



'' ^^<-^'' (31) 



15(2)i87r5a3d-10' 

b. The term with two derivatives with respect to the components of p may be calculated 
analogously. One finds 
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such that the pole part is 
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d. Similarly, the term without derivatives computed at zero external momentum furnishes the 
pole part 

• ^ ^ ^ (-56^ + 4a2m2). (36) 



Collecting the above results we can write 



II. To complete our computation of the counterterms, we need to consider the three point 
function. At one-loop a direct calculation gives that 

1 A"^ 1 

r^^) = -iX + i ——- — \- finite terms . (38) 

We are now in a position to determine the renormalization group flows of the parameters 
of the model. In fact, the A^ point vertex functions of the model satisfy a 't Hooft- Weinberg 
renormalization group equation 
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r^^\p,m'^,b'^,a'^,X,f,) = 0, (39) 



where the renormalization scale /i was introduced by replacing A — t- /i'^"A, where e = 10 — d, 
so that the new coupling constant A is dimensionless. After removing the pure pole part, the 
renormalization group parameters may be then obtained by replacing into (|39|) the renormalized 
vertex functions 

r(2) = i[pl - 62p2 _ a2(p2)2 _ ^2 _ ^2(-p-j^-^g^ _ i^ /^Residuei)] (40) 

and 

r(3) = -iX + iA^ (Finite2 - In /xResidue2) , (41) 

where Finitei^2 are the finite parts and Residuei_2 the residues at the poles of the corresponding 
vertex functions. We find 
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which agree with [5|] and shows that the model is asymptotically free. Furthermore, the effective 
mass increases or decreases with fj, accordingly b^ is greater or less than 2/3?n,^a^. By introducing 
a logarithmic scale, t = ln(;u/;Uo), where /io is a reference scale where the parameters in the 
Lagrangian have been defined, we may evaluate the flows of the effective parameters as follows, 

r2 
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These equations imply that 



with a(0) = a (43) 



with A(0) = A. (44) 
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whose general solution is 



m = - — ^^-^ , (46) 

(19t + 6ci)T9 
where ci and C2 are constants determined by the initial conditions. Using this result, we integrate 
the equation for a giving 

a^(t) = ^(19t + 6ci)^ (47) 

and also 
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Notice that a decreases with t vanishing at the critical value tjR = —^wti- The same happens 
with the effective parameter b which is given by 

f (t) = ^ ^ (19t + 6ci ) ^ . (49) 

(6ci)i9 

However, at the same time the effective coupling constant increases tending to infinity as t ap- 
proaches tin. Our perturbative methods are not applicable insofar A is not small and no conclusion 
can be made concerning the restoration of Lorentz symmetry at small momenta. 



1. Unitarity in the anisotropic situation. 

Differently from what happens in Lorentz invariant theories but with higher space and time 
derivatives, models where only higher spatial derivatives are introduced preserve unitarity (in this 
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respect see also [9|). Let us exemplify this by analyzing the Cutkosky rules [10|] for the one- 
loop contribution to the two point function of the model (|26p . Those rules demand that the 
diagrammatic relation of Fig. [2] be obeyed. The analytic expression for the graph on the left hand 
side is 
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where the only restriction on w^ is that it depends only on the spatial part of the momentum. For 
the model ([26]) . ujf, = 6^k^ + a^(k^)^ + m?. By integrating over k^ and using the identity 

' ^" (51) 



where P denote the Cauchy principal value, we get 

On the other hand, the expression for the cut diagram on the right of Fig. [2] is obtained by the 
following replacement of the free propagator 
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and by noticing that the vertex factor on the left of the cut is the complex conjugate of the one on 
the right of it. Therefore, the graph on the right of Figl2] gives 



which after integrating over k^ gives the same result as in Eq. ([52]) . 
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FIG. 2: Cutkosky rule at one-loop order. 
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IV. THE ^^ MODEL 

The if^ model with z = 2, specified by the Lagrangian density 

turns out to be renorniahzable in 6 spatial dimensions, the degree of superficial divergence being 
given by 

d{G) =8-2N. (56) 

If a BPHZ like scheme is adopted, we find that the Green functions satisfy the equation of 
motion 
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where X = Y\i<f{xi). This expression may be derived by noting that in momentum space the 
operator applied on ip, in the left hand side of the above equation, is equal to —i times the inverse 
of the free field propagator 

To fix the renormalization group parameters, we compute the radiative corrections to the pa- 
rameters of the model as follows. 

The lowest order correction to the two point function comes from the tadpole graph in Fig. [3] 
whose analytic expression is 

A f dko d'^k 1 



Since the above integral does not depend on the external momentum, with the help of (|A6P its 
divergent part is easily obtained, i.e., 

.(-362+4a2m2) A 1 
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The one-loop contribution to the coupling constant renormalization it is also straightforwardly 
calculated. In fact, we just need to calculate the pole part of 
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FIG. 3: Two point function of order A. 
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which, again with the help of ()A6p . is found to be 
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To find the lowest order corrections to the a and b parameters is a much more difficult task since 
it involves the calculation of a two-loop graph. Actually, only for 6 = and m = closed analytic 
expressions exists ll[. In that case one needs to calculate 

i\^ f dko d'^k dqo d'^q I 1 1 
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3! J 2tt (27r)'^ 2tt (27r)<^ k^ - a^{k^y ql - a2(q2)2 {p^ - ko - q^f - a?{{p - k - q)2) 
whose pole part turns out to be 
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yielding the following beta function 
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Similarly, the computation of the pole part of the four point vertex function furnishes 
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so that the model is infrared stable. This result will be used in the next section as the renormal- 
ization for the Yukawa model demands the inclusion of the ^p^ self-interaction. 

Concerning to the solution of these equations, there is a particular configuration of the initial 
values of the parameters a and A, such that the system of differential equations above exhibits 
a simple analytical solution. In fact, denoting the constants appearing in above equations by 
A = j^ (II -|- 31n I) and B = ^^^, we can get from ([65]) and (f66|) . the following equation for the 
effective coupling constant 

3A rd\\^ „^/5A\^ ^2 (d'^X\ 



'5A\'_ -/5A\' ^2 
2^2 vatj [dtj + \dt^J 



0, 



(67) 
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3^t 



which admit the solution 

A(t) = \{0)et^\ (68) 

From (|65p . we obtain 



a{t) 



2 



9A 

Ib2 



X \0){e^'-l)+7f{0) 



(69) 



l2-r2. 



But this win be a solution of the system only under the condition a^(0) = j^X"{0), yielding 



1 



a(0 = (I^a'(O)J e^*. (70) 

This result show us that the Lorentz symmetry breaking parameter a goes to zero as t — t- — oo, 
i.e., in the infrared region. On dimensional grounds, we see that these solutions would not suffer 
modifications even in the presence of the parameters b and m, such as in lagrangian (I55p . So, we 
concludes that in this case the Lorentz symmetry can be approximately recovered in a sufficient 
low scale of the energy. 

V. A RENORMALIZABLE YUKAWA MODEL WITH z = 2 

The analysis of the previous sections indicates that anisotropic models that are infrared stables 
are the best candidates to show Lorentz symmetry restoration at small energies. To further inves- 
tigate this possibility we consider now a model of boson and fermion fields interacting through the 
Lagrangian density 

1 6^ rn^ c? 

C- = -doifdoip - -^dupdup - —Lf^ - -^ipA'^ip 

+ V5(i7°5o + ib^j'di + a^A - M)^; + igi^T^^jip - ^/, (71) 

where T^ is the chiral matrix. In momentum space, the free propagators derived from the above 
expression are 

1 

^^^^^ = A;2-62k2_„2(k2)2_^2' (^2) 

for the bosonic field ip and 

"^^^^^ = 70^0- 6^7 ■k-a^k2-M' ^^^^ 

for the fermion field. The interaction between the basic fields is given by the non-quadratic part 
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FIG. 4: Interacting vertices for the Yukawa model. 

of (17ip so that the corresponding Feynman rules are as depicted in Fig 21 With the anisotropic 
scaling exponent z = 2, the model is renormalizable in d = 6 spatial dimensions and one time 
dimension. Indeed, in this situation the degree of superficial divergence for a generic graph G is 

d{G) = 8 - 2Nb - 3Nf, (74) 

where Nb and Np are the number of external boson and fermion lines of G. As the spacetime 
dimension is odd, we will work with gamma matrices adequate for eight dimensions as to be able to 
define the chiral matrix as (notice that this assumption is in accord with the fact that the effective 
dimension of the underlying spacetime is eight) 

7 
At=0 

The above choice for the interaction term prevents the induction of counterterms proportional to 
ip^, if and ■i/'^c^ and so is simpler than the other possible renormalizable interaction tptpcp. 
As usual, we define renormalized quantities through the replacements 



^ ^ 4/V, (76) 

^ ^ 4^V, (77) 



9 ^ -^^S, (78) 

A ^ I^A (79) 



^r 



and also the Lorentz breaking parameters. 



{bl,b^,al,a^} -^ {{Zb^/Z^)bl,{Zb^/Z^)b^,{Za^/Z^)al,{Za^/Z^)a^}. (80) 

For the computation of the renormalization constants, we shall use the dimensional reduction 
scheme I4] in which all algebraic simplifications are done with the Dirac matrices as introduced 
above and, afterwards, the integrals are promoted to d = 6 — e spatial dimensions. The ambiguities 
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that may be present in such procedure [13[] only manifest themselves in higher order and do not 
affect our one-loop computations. Up to one-loop order, the relevant graphs are shown in FigslMSl 

b 
a 




FIG. 5: Lowest order contributions to the self-energy of the scalar field. 



We begin by considering the order g'^ radiative correction to the self-energy of the field 99, Fig. \5^, 
whose unsubtracted analytic form is 



^M = 9' I 



dko d'^k 

'2tt {2-kY 



Tr 



f 



-i^k^ - 6^7 ■ k - a^k2 - M 



X r^ 



(81) 



^0(^0 + pO) _ 5^^ . (k + p) - a^,(k + p)2 - M_ 

where, as we are considering the eight dimensional representation, the gamma matrices are 2 
dimensional and obey Tr7'^ = and Tr7'^7'^ = 2'^g^'^. 

As before, for generic momenta the computation of S^(j)) is cumbersome and probably unfea- 
sible. However, to extract its pole part we may proceed as in section IIIII After discarding some 
finite contributions, the term of the fourth order derivative with respect to the spatial momenta 
gives 



2r,4: PiPjPkPl 

^ 4! 



dko d'^k 

'2tt {2ttY 






d^ 



1 ^al\^ 



d^ 



k2 



[den] dkidkjdkkdki [den] [den] dkidkjdkkdh [den] 



(82) 



where, for simplicity, we have defined den = k^ — hV^i? — (a^k^ + M)^. Thus, we have to compute 



integrals of the form 



J{x,y,z) 



dkn d'^k 



Kq|K| 



(83) 



2tt {2'kY [kl - 62 k2 - (a^k2 + MfY ' 
which, up to some finite terms, can be read from the appendix El by sequentially making the 



replacements rm? — )■ M2, iP' ^ h\ + 2Ma^ and a? -^ a\- The pole part at (i = 6 gives then 



(2)77r3 a^d-Q 



[p'?. 



(84) 
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The other terms in the Taylor expansion may be calculated similarly. We find that the pole part 
of the self-energy of the (p field is 



(364 + 12a^62 ^ + 8a2 m2) ^ ^ ^ 



^4> 
647r3 



< 32vr3a2 



Po + 



M 1 



P^ + 



1 



2\2 



r:;3(p^) 



5 



1 



a^ d — 6 



(85) 



> 3(2)7^3 

The tadpole graph has already been considered in our study of the 99^ model and we simply quote 
that the result is given by Eq. ([60]) . 

Let us now consider the one-loop contribution to the self-energy of the fermion field which 
corresponds to the graph in Fig. [6l The analytic expression is 

1 



^i>{p) 



2 f dko (fk 5 

'^ J 2tt {2ttY l^k^ - 6^7 • k - a^k2 - M ^ 

1 

' {PO - koY - blip - k)2 - a2 [(p - k)2]2 - ^2 

2 / dko d'^k -7°/c° + 6^7 • k + a^k2 + M 



-g 



27r {2ttY kl - 52 k2 - (a^k2 + Mf 
1 



(po - koY - 62 (p - k)2 - a2 [(p - k)2]2 - m2 ■ 
Prom this we get the pole terms: 



(86) 



FIG. 6: Two point fermionic function. 



a. Term with two derivatives with respect to the spatial part of the external momentum: 

-i a^(3ay + a^) g2p2 
3(2)'' 7r3 a^{a^ + a^Y d - 6' 

h. Term with one derivative with respect to pq: 



9 



-POT 



(2) ''vr^ a^[a.^ + a^Y d-6 
c. Term with one derivative with respect to the spatial part of the external momentum: 



(87) 



(88) 



(2a<^ + a^)b^ g^ 



3(2)'^ 7r3 a^a^{a^ + a^Y d-6 



P-7- 



d. Term without derivatives: 



1 ^2 [a2(2a^ + a^)62 +a2(a^ + 2a^)62 +2a2a2M] 1 



2^3 ^3^2 



(a^ + a<^)" 



d-6 



(89) 



(90) 
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FIG. 7: Three point vertex function. 



The one-loop contribution to the three point vertex function, Fig. [71 furnishes 



3^5 



:9'^ 



whereas for the four point vertex function of the bosonic field shown in Fig. [8] we obtain 



i , 6 



1 



IGvT'^ a' 



-A^ 



1 






The vertex functions of this model satisfy the renornialization group equation 



(91) 



(92) 



\ 



>' 






\, 



FIG. 8: Four point scalar functions. 



^'djl^^'^^^^^ ^3-Q^ + ^^dX ~ ^""^^ ~ ^^^^ 



where, we have defined the differential operators 



r(^)(p,m2,6,a,A,/i) = 0, (93) 



d d d 

and 



(94) 



^^='^M^+/?6,^+/3a,^ 



"(/j 



(95) 



20 



The parameter fi was introduced through the following replacements of the coupling constants 
g — )■ gfi^ and A — > A/i*^. 

By substituting the one-loop expansions of the renormalized vertex functions we then get the 
various beta functions (see the appendix [B] for details), with the arguments indicating the order of 
the coupling constant in which they have been calculated: 



/?6^,(5^ 



1 {hl + 2a^M) 

—Ji 9 > 



327r3 



V 



(96) 



f^aiig' 



1 (4^ + 4) 2 



1287r3 



(97) 



2^ ^ fey. (fly) + 2a^) 2 
'''^^ 1927r3 a<^a^(a<^ + a^)2 



l^b^ia 



(98) 



P-A9 ) 192^3 a^(«^ + «^)3 5 ' 



(99) 



and 



f^gig') 



2567r3 



2a3 + ia^a^ + 3a^a^ + 2a^ 



/3a 



5127r3 a3 



3 ff^ 1 

+ 



Xg' 






167r3 



(100) 



(101) 



As the Eq. (|96p involves the mass of the fermion field, we need also the corresponding renormal- 
ization group function: 



[2a^alhl + alhl + albl + 2ala^{bl + 2a.^M)] ^ 

ouyg ) = nr^_?!.,:^, t. — r~~~:^ g 



(102) 



2bQ-K^d^a^{a^ + a^y 

The equations that govern the evolution of the effective Lorentz symmetry breaking parameters 
may be separated in two sets. In the first one are the equations for a^p, a^ and g, which do not 
depend on the remaining parameters. In the second set are the equations for the other parameters 
b^ and b^, those need the input of the former set to be evaluated. Unfortunately these equations 
do not seem to have simple analytic solutions so that we use numerical methods to investigate 
their properties. Before proceeding, we need to stress what we mean by restoration of the Lorentz 



21 



symmetry in this model. Obviously, we can never take the symmetry breaking parameters a^ and 
a^p equal to zero, since we would end up with a nonrenormalizable theory. In this sense, we can not 
expect an exact Lorentz symmetry, which would correspond a,^ = a^ = and b^ = b^; these are 
not even special points (for example, fixed points) of the above renormalization group equations. 
However, we can consider the possibility of an approximate Lorentz symmetry arising in some 
specific low energy limit, depending on a fine-tuning of the parameters involved. More precisely, 
we may find a region where the parameters o^ and a^ are sufficiently small, and, furthermore, the 
parameters b^p and b^ are such that b^p k, b^. Of course, this is not the ideal situation, but it could 
furnishes a positive view concerning the Lorentz symmetry restoration in anisotropic field theories. 
To get a better insight on this possibility, we performed a numerical study, as described bellow. 

Our results, obtained through the use of the numerical package of the Mathematica (and also 
checked with the Runge-Kutta 4th order method), show that the fermionic parameters change 
much more slowly than the bosonic ones. Actually, the changes of the tangents to the curves of 
some parameters are so small producing the impression that they are straight lines. At higher 
momenta, the coupling constant 'g, a^ and b^ increase steeply for t ~ 10^, indicating the existence 
of a singularity similar to the Landau pole found in many not asymptotically free field theories, see 
Figs. [9lll0landlll[ By contrast, a^ and b^ increases slowly as shown in Fig. [QJandfTOl For negative 
t, corresponding to the small momenta region, the general pattern is that all parameters decrease, 
as remarked before, the parameters associated to the ^ field do that in a more slow rate than those 
associated to the ip field (see Figs. [T2] and [T3|) . The graphical figures were drawn by taking the 
initial value ^(0) = 10"'^ for the coupling constant and the initial values of the other parameters 
all equal to one. Variations of these initial values do not qualitatively change the behavior of the 
effective parameters. It should be noticed that below the lower ends of the curves the data are 
not reliable as the modulus of t is very large whereas, at the same time, the effective parameters 
become very small possibly generating large numerical errors. 

Thus, given a low energy region, in order to have an approximate Lorentz symmetric situation, 
one should modify the initial configuration of the parameters o^ and a^, so that they will attain 
the desired range of values in the specified region. For example, as the a's monotonically decrease 
when t is negative, this can be easily done by choosing the initial values of a^p and a^ to be equal 
to the maximum values they should have in the region of interest. Moreover, we need to adjust 
the initial configuration values of the b^p and b^ parameters, so that 6<^ ~ b^ in the same region. 
Notice that, these two requirements can be satisfied due to the decoupling of the equations for a's 
and 6's, as mentioned before. 
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FIG. 9: The ultraviolet behaviors of the parameters a^p and a^. 
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FIG. 10: The ultraviolet behaviors of the parameters b^ and b^. 
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FIG. 11: The general behavior of the coupling constant g. The singularity in the UV region is similar to 
the Landau pole. 

VI. CONCLUSIONS 

In this work we have examined some aspects of anisotropic field theories as renormalization 
properties and Lorentz symmetry restoration. In particular, we performed a renormalization group 
analysis of various models aiming to understand the infrared behavior of their effective parame- 
ters. Initially, we considered models without higher derivatives but with fields with different light 
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FIG. 12: The infrared behaviors of the parameters a^ and a^. 
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FIG. 13: The infrared behaviors of the parameters 6^ and b^. 

velocities; then in two specific instances we verified that the breaking of Lorentz invariance is soft, 
i.e., the Lorentz symmetry is recovered at low energies. 

In the general situation where higher spatial derivatives are present, unless for super- 
renormalizable models, the idea that Lorentz symmetry is restored at low energies requires that the 
interactions be infrared stable and some process of dimensional reduction, to cope with eventual 
divergences appearing whenever the higher derivative terms are eliminated. It should be pointed 
out that, as may be straightforwardly derived from Eq. (j22p . the spatial dimensions in which the 
anisotropic, d, and Lorentz symmetric, di^ versions of a given model without derivative couplings 
are renormalizable are related hy d = zdi where z is the critical exponent which characterizes the 
anisotropy. 

In the case of the 99'^ model, although the breaking parameters decrease with the energy, because 
of the asymptotic freedom of the model, the effective coupling increases and no conclusion can be 
achieved on the restoration of Lorentz symmetry within this perturbative approach. 

Concerning the behavior of the parameters under the renormalization group a more favorable 
situation occurs for the ip^ and Yukawa models which are infrared stable. For the 93^ model, we 
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found a special analytical solution in which the effective parameters increases monotonically from 
zero in the infrared to very large values at high momenta. For the Yukawa model the system of 
equations that govern the changes of the effective parameters are intricate enough and no closed 
analytical expression seems to be feasible. Our numerical solution revealed that for high momenta a 
singularity, like the Landau pole, is present but for small momenta all Lorentz breaking parameters 
decrease. As discussed in the text, our results indicate that we can find specific low values of the 
energy scale where the approximate Lorentz symmetry may be achieved (in the sense of small o's 
and bip ~ b^), what nevertheless requires a fine-tuning of the initial values the parameters. 

In this study we have restricted ourselves to models with z = 2. As we have seen, the inves- 
tigation of the Lorentz symmetry restoration involves the determination of the infrared behavior 
of the Lorentz breaking parameters and thus for models with higher values of z, which have more 
parameters, the situation becomes more complex. Another lateral remark concern models with 
gauge symmetry; in this case, to keep the symmetry, covariant derivatives have to be used so that 
in the perturbative approach new interaction terms must to taken into account. For a discussion of 
the symmetries using the BPHZ approach on anisotropic models, see the sequel [ij] of this work. 

We hope that the ideas and methods we presented may be useful also in condensed matter 
physics, in the contexts of quantum phase transitions and Lifshitz models, i.e., in situations where 
the kind of anisotropy here considered may be a natural concept jlSl. Il6|. 
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Appendix A: Integrals 

In this Appendix we will compute the pole part of some integrals that are relevant in our 
calculations. We begin by considering the following integral 

J{x,y,z)-J 2vr (27r)'^[/t2-62k2_a2(k2)2-m2 + ie]-' ^ ' 

where x, y and z are such that the integral is at most quartically divergent when the dimensional 
regularization is removed. 
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We use the representation 
J{x,y,z) 



i^T{z)J 27r(27r)^^°'^' io 

\z) io y 27r ° y (2^)'^ ' ' 



i^[_62k2_„2(,^2-)2] 



i^r 



Now, 






^fcje^^^o = J_(l + (_i)^)(_i^)-^^r 



{l±£i^,X + 1, 



and, denoting by Jl^ = 27r°'' ^/r(ci/2) the volume of the (i dimensional unit sphere, 



(A2) 



(A3) 



d'^k 



IJ^iy n[-62k2_„2(k2)2] ^ ^d (^Q 7) 

(27r)'^' ' (27r)'^ 4 

Nl/2 



2^N-(d+2/)/4 



■ ,(i + y, /^c? + y 1 «7^^ 

J- ( — : — J 1-^1 



4 '2'4a2 



(n)'/ ,2T..2 + t^ + y), ^ (2 + d + y) 3 ^76' 



-6^r( 



Fi 



, (A4) 



\a\ • 4 ' ^'^ \ 4 '2' 4a2 

where iFi is the confluent hyper geometric function. The divergence that appears in ()A2p when 
the dimensional regularization is removed is due to the behavior of integrand for small 7. Thus, to 
obtain its pole part we may use the approximations, 



iFi 



d + y 1 ijb'^ 



i-t'i 



4 '2'4a2 
2 + d + y) 3 i-fb^ 



1 + i 



b\d + y) 



8a2 

1 + 0(7). 



7 + 0(7'), 



(A5) 



4 '2'4a2 

Notice that in the last expansion we need to consider just the zeroth order term because of the 
7^/2 additional factor in ()A4p . All powers of 7 greater than those that we have considered produce 
finite results. It is now straightforward to perform the remaining integrals and we obtain 



J{x,y,z) 



1 {l + {-ir)^fx + l\ a 



i^r(z) 



,-(d+y+4) 



(4vr).+ir(|) 



a^T 



d + y 



+ iV 



X r 



d + y\b'^ 



:(d + y)(_l)|(2-<i-J/)e^(,^2)l(-2+d+2.+y-4.) 



-2 - d - 2x - y + 4z 



2-d-2x-y + Az 



iaH^T 



4^ /d + y + 2 



X (im2)3('^+2-+2^-4-)r 



4 J V 4 

d — 2x — y + 4z 



(-1) 



U-d-Ux-y) i2TTX 



(A6) 



Another typical integral that occurs in our calculation involves the product of propagators with 
different parameters as in 



dko d'^k 



27r (27r)rf(A;2-62k2-a2(k2)2 



miY^ {ki - bik^ - 4{k^y 



mlY-i ' 



(A7) 
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In this case, we first join the two denominators by using Feynman trick 

A^iB^^ T{zi)r{z2) Jo [Ax + B{l-x)Y-^+^2 ^ ' 

and then apply the formula ()A6p . 



Appendix B: Renormalization group parameters 

Here we will present some details of the derivation of renormalization group parameters listed 
in the text. First notice that the boson and fermion two point vertex functions are given by 

^\^ \P) = ^bo ~ ^v^P^ ~ o^(P^)^ ~ "^^ ~ 5^(Finitei — ln//Residuei) 

— iA(Finite2 — ln/iResidue2)], (Bl) 

^f{p) = ^bV - 6^7 • P - avP^ - ^ + ^^(Finites - In/iResidueg)], (B2) 

where the residues have the form 

Residuei = ^i + ^2Po + ^sP^ + ^4(P^)^ and Residue2 = M, (B3) 

Residues = -Bi + S27 V + -^sT • P + ^4 P^ (B4) 



and the A's and i?'s can be read directly from (j85p and (J87ti90p . Similarly, the three point fermion- 
boson vertex function and four point of the boson field have the expressions 

r(3) = -5r^-53(Finite4-ln^Residue4), (B5) 

r^"^^ = —i\ + (7^(Finite5 — In/xResidues) — A^(Finite6 — In/xResiduee), (B6) 

where Residue j, i = 4, 5, 6, are given by 

Residue4 = C,= 1,3,3 ,^,2 | ,2 ,/ ^ I^-idue^ = D, = ^^^^, (B7) 

Residues = 5i = z^-^. (B8) 



By replacing these expressions into (|93p and equating to zero the coefficient of each power of the 
coupling constant we get 

/3a (5') = iD^9\ /3a(A2) = -ib^X^ (B9) 

and 

/3a(A5') = 4A7^(5'), (BIO) 
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so that 



/3a = ^-- - — ^ + ^-^ 



(Bll) 



We have also 






(B12) 



/?a2(5' 



1 (4a2 + 0,2 



(«^^^ + ^4).^-I^ «■ 



3 y ' 



V' 



(B13) 



''' 1927r3 0(^0^(0;^ + a^)2 



(B14) 



and 



/3,,(<7 ) - -{a^B, + i?4)5 - i92,3,^(„^ + „^)3 5 



(B15) 



/3.(5' 



2 ^ ^ 1 y y 2567r3 



2a3 + 40^0^ + 3a<^a^ + 2a^ 



(B16) 
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